Matter rogue wave in Bose-Einstein condensates with attractive atomic interaction 
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We investigate the matter rogue wave in Bose-Einstein Condensates with attractive interatomic 
interaction analytically and numerically. Our results show that the formation of rogue wave is mainly 
due to the accumulation of energy and atoms toward to its central part; Rogue wave is unstable 
and the decay rate of the atomic number can be effectively controlled by modulating the trapping 
frequency of external potential. The numerical simulation demonstrate that even a small periodic 
perturbation with small modulation frequency can induce the generation of a near-ideal matter 
rogue wave. We also give an experimental protocol to observe this phenomenon in Bose-Einstein 
Condensates. 
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I. INTRODUCTION 

The dynamics of Bose-Einstein condensates (BECs) 
at ultralow temperature are described well by Gross- 
Pitaevskii (GP) equation [l|, in which the nonlinear- 
ity is arose from interatomic interactions characterized 
by the s-wave scattering length. Recently experiments 
have demonstrated that tuning of the effective scattering 
length, including a possibility to change its sign, can be 
achieved by using the so-called Feshbach resonance tech- 
nique [2]. In particular, the experimental realization of 
BECs in dilute quantum gases have opened the floodgate 
in the field of atom optics and condensed matter physics 
At the same time, the collective excitation of matter 
waves in BECs has also drawn a great deal of interest to 
explore the dynamics of BECs deeply from both experi- 
mental and theoretical perspectives, such as matter wave 
solitons 0-01, periodic waves jTo[ , shock waves [11], vor- 
tex [l2j and necklaces [l3j|. However, to our knowledge 
less attention have been paid to the matter rogue wave 
which is a fundamental and novel nonlinear excitation in 
BECs. 

Rogue waves from ocean are that their heights, from 
crest to trough, are more than about twice the signifi- 
cant wave height [14]. They appear without any warning 
and disappear without the slightest trace. Owing to se- 
vere environment and high risk in ocean, the systematic 
study of rogue waves become so difficult that the nec- 
essary conditions and physical mechanism of their gen- 
eration are not sufficiently well understood. Recently, 
theoretical studies have shown that the rogue wave phe- 
nomenon can be explained well by nonlinear theories 
[HI [HI , and the various possible formative mechanisms 
have been discussed, such as the modulation instability 
in one dimension [I?], [H| , nonlinear spectral instability 
[l9| and in two-dimensional crossings j2Cj. Furthermore, 
the rogue waves phenomenon have been observed exper- 
imentally in variety of physical systems including optical 
fibers [U, [22| , arrays of optical waveguide [23[ and cap- 



illary waves [24[ . 

As a nonlinear physical system with similar nonlinear 
characteristics, BECs can support the interesting rogue 
waves and allow us to understand deeply the nature and 
the dynamics of rogue waves in laboratory conditions. 
The management of Feshbach resonance for nonlinearity 
and a tunable atomic trapping potential also provide us 
with a powerful tool for manipulating rogue wave. This 
enable BECs to have more advantages for investigating 
rogue wave than other physical media. In this paper, we 
mainly investigate the matter rogue wave of Peregrine 
type with the emphasis on its formative mechanism in 
BECs. We firstly obtain the exact rogue wave solution of 
the GP equation with time-dependent attractive atomic 
interaction in an expulsive parabolic potential. By an- 
alyzing the atomic number density distribution in the 
rogue wave against the background, the formative mech- 
anism of matter rogue wave can be clarified that the ac- 
cumulation of energy and atoms toward to its central 
part. Rogue wave can not keeps dynamic stability and 
the decay rate of the atomic number can be effectively 
controlled by modulating the trapping frequency of ex- 
ternal potential. Moreover, we use the breather evolution 
in a regime approaching the excitation of matter rogue 
wave as approximation to simulate the creation of mat- 
ter rogue wave numerically, which indicates that a small 
periodic perturbation with small modulation frequency 
can excite a near-ideal matter rogue wave. Finally, we 
also give a practical and effective experimental protocol 
to observe this interesting phenomenon in future BECs 
experiments. 



II. MATTER ROGUE WAVE SOLUTION 

Under the mean-field level, the evolution of macro- 
scopic wave function of BECs obey the 3D GP equation 
[3| . For a cigar-shaped condensate at a relatively low den- 
sity, when the energy of two body interactions is much 
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FIG. 1. (Color online) The asymptotic processes from Eq. j2]) 
to Eq. J3} in the limitation of (A s ,k s ) — ^ (2A c ,/c c ). As the 
bright soliton amplitude A s and wave number k s approaching 
(2A c ,k c ), the spatio-temporal separation between adjacent 
peaks gradually increases shown in Fig. l(a-f), where the 
parameters are (a) A s = 2.4, k s = 1.2. (b) A s = 2.2, k s = 
1.2. (c) A s = 2.01, k 8 = 1.01. (d) A s = 1.6, k s = 0.8. (e) 
As = 1.8, k a = 0.8. (f) As = 1.99, fc s = 0.99. Fig. 1(g) 
represents the exact matter rogue wave solution Eq. (3} with 
A s = 2, k s = 1. Other parameters are A — 0.01, to — 10, k c — 
A c = 1. The red and green lines represent the sloping lines 
Vq and V«, respectively. 



less than the kinetic energy in the transverse direction, 
the system can become effectively quasi-one-dimensional 
regime with time-dependent attractive interaction in an 
expulsive potential |25j . 



1 d 2 ^ 
2~dx^ 



+ a(t)|^| 2 ^+^A 2 x 2 ^ = 0, (1) 



where the aspect ratio reads A = \ujo\ /w± <C 1, coordi- 
nate x and time t are measured in units a± and 1/uj± 
with a± = ^h/muj±_ (m is the atomic mass) and ao = 
y/h/mujo the linear oscillator lengths in the transverse 
and cigar- axis directions, respectively. uo±_ and uoo are cor- 
responding harmonic oscillator frequencies. The nonlin- 
ear coefficient a(t) is defined as a(t) = \a 8 (i)\/aB, where 
a s (t) is so-called s-wave scattering length. Correspond- 
ing to the real BECs experiment [9[ , the bright soliton can 
be created for 7 Li by tuning the scattering length contin- 
uously in an expulsive potential with uj± = 2tt x 700Hz 
and ujo = 2m x 7 Hz. So we can choose the nonlinear co- 
efficient in the form of a (t) = exp[A(t — to)] manipulated 
by Feshbach resonance technique [26|, where to repre- 
sents an arbitrary real constant determining the initial 
scattering length \a s (t = 0)| = aBe~ xt °. 

To obtain the exact solution of Eq. (pp), we introduce 
the transformation ip = q(X, T) exp[A(t — to)/2 — i\x 2 /2] 
with the coordinate transformations X = e x ^~ to ^x and 
T = [e 2 ^-^ - 1]/(2A). Then Eq. (TJ) can reduce to 
the standard nonlinear Schrodinger equation, and the so- 
lution of Eq. (pQ) constructed on continuous wave few) 
background i(j cw — A c e lip can be obtained as follows [27j, 



. y cosh + cos a A r\ sinh 6 + 5 sin a 
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where (9 = M/X - [A S M^ + (fc c + fc s ) M 7 ]T/2, a = 
M^X - [(fc c + fc fl )Mii - A a M/]T/2, <p = <p c - Xx 2 /2 - 
i\(t -t )/2,(p c = k c X + (A 2 - /c 2 /2)T and Mr + zMj = 
V^(fe c - fc s - z'A,) 2 + 4^2. The subscripts R and / de- 
note the real and imaginary part of M, respectively. 

In general, some excited state solutions, such as cw 
wave and bright soliton on the background of ground 
state, can be recovered successfully from Eq. (|2j). Firstly, 
when the cw background amplitude vanishes (A c = 0), 
Eq. ([2]) reduces to the bright soliton solution ip so i = 
A s e liPs sech9 s with varying amplitude A s e x ^~ to ^ and 
group volecity v s = k s cosh[A(£ — to)] [28], where S = 
A s (X-k s T) and <p 8 = k s X + (A 2 - k 2 s )T/2, and k s repre- 
sents the wave number of bright soliton. Secondly, when 
the initial amplitude of bright soliton vanishes (A s = 0), 
Eq. ([2|) reduces to cw background solution i(j cw = A c e lip 
with varying group velocity v c = k c cosh[A(t — to)], where 
k c is the wave number of cw background. So the exact 
solution Eq. ([2j) represents a bright soliton embedded 
in a cw background field. It should be pointed that the 
parameters A C1 A S1 k c and k s are derived from the math- 
ematical construction of Eq. (2), the numerical values of 
these parameters can be chosen freely, so we assume that 
these parameters are real constants without loss of gen- 
erality. Furthermore, the dynamical evolution of solution 
Eq. ([2]) are shown in Fig. l(a-f). From Fig. 1 we observe 
that the solution in Eq. ([2j) commonly exhibits a breather 
characteristic and a time periodic modulation of the soli- 
ton amplitude, which can be regarded as the results of 
the interaction between the localized process of cw back- 
ground along the slope direction V a and the periodiza- 
tion process of bright soliton along the slope direction V#, 
i.e., the bright soliton undergo periodic energy and atoms 
exchange with cw background, where Vq and V a repre- 
sent the lines MjX - [A S M R + (k c + k s ) M/]T/2 = and 
M R X - [(k c + k s )M R - A s Mi]Tj2 = on the space-time 
plane as shown in Fig. 1, respectively. 

Through above analysis we can see that when k c = k Sl 
the bright soliton must move with the same group ve- 
locity as that of cw background which implies that the 
bright soliton and cw background are in "resonant state" 
in space. The critical point A 2 = AA 2 forms the dividing 
line between modulation instability process (A 2 > A 2 /4) 
and the periodization process of bright soliton (A 2 < 
A 2 /A) under the resonant condition k c = k s [27[, which 
means that the relative initial intensity between the cw 
background and bright soliton determine the different 
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FIG. 2. (Color online) The plots for the evolutions of a\ and 
CX2 as the function of A s and k s during the limit processes, 
where k c = 1 and A c = 1. The two angles ai and «2 as the 
function of A s and k s are not continuous at the limit point. 



physical behaviors of the solution Eq. (j2]). Especially, 
with the limit conditions of k c = k s and A 2 S = 4A^ we 
have 



4 + iSA 2 r T 



1 + 4A*T 2 +4A 2 (X-k c TY 



A c e v 



(3) 



which represents a localized matter wave with the max- 
imal amplitude Ap — 3A C in BECs, and its dynamical 
evolution is shown in Fig. 1(g). Interestingly, the exact 
solution Eq. ([3]) displays the typical rogue wave charac- 
teristics of Peregrine type that a localized breather char- 
acteristic with only a single hump both in space and time, 
which indicates that the localized wave is captured com- 
pletely at x = and t = to by the cw background 
So far, such solution Eq. (|3j) has been conjectured to be 
a prototype of oceanic rogue waves. 



III. DYNAMICS OF MATTER ROGUE WAVE 

In order to better clarify the formative mechanism of 
rogue wave solution in Eq. (j3j), we firstly investigate the 
asymptotic processes of Eq. ([2j) to Eq. ((3]) in the limit 
processes (A s ,k s ) — »• (2A c ,k c ) by fixing the numerical 
values of cw background amplitude A c and wave num- 
ber k c . From Fig. l(a-f), we can observe clearly that the 
spatio-temporal separation between adjacent peaks grad- 
ually increases as the bright soliton amplitude A s and 
wave number k s approaching (2^4 C , fc c ), which leads to a 
greater spatio-temporal localization in Eq. ((2]). Further- 
more, the parameters tanai = 2/(k c + k s + A s Mr/Mi) 
and tan«2 = 2/(k c + k s — A s Mi/Mr) shown in Fig. 1(c) 
represent the slope of the lines Vq and V a at x = and 
t = to, respectively. When the values of A s and k s ap- 
proaches to the critical point, Vq and V a gradually turn 
to a relative fixed direction associated with the maxi- 
mal spatio-temporal localization in Eq. (j2j). However, 
the two angles a\ and ot2 as the function of A s and k s 
are not continuous at the limit point shown in Fig. 2, 
i.e., ai and ot2 do not exist limitation, which play the 
important role to describe that the solution Eq. (j2j) is 
localized along Vq and V a at the limit point. Especially, 
rogue wave solution in Eq. (j3|) can be considered as a 




FIG. 3. (Color online) (a) The atomic number density dis- 
tribution in matter rogue wave at fixed time point, (b) The 
atomic exchange between matter rogue wave and background. 
The parameters are A c = k c = 1 and to — 10. 



transition state between the modulation instability pro- 
cess (A s — » 2A~) and the periodization process of the 
bright soliton (A s — >> 2A+) under the resonant condition 
kc — k s . 

As shown in the following, the formation of rogue 
wave can be clarified by the atomic number density dis- 
tribution against the background defined as p(x,t) = 
\t/j R (x,t)\ 2 - \ip R (x = ±oo,t)| 2 . With Eq. ©, we have 



SA 2 C + 32A 6 C T 2 - 32A*(X - k c Tf 
[1 + 4A 4 C T 2 + 4A 2 (X - k c T) 2 } 2 



3 A(t-t ) 



(4) 

and the time-independent integral f_™ p(x, i)dx — 0. 
From the condition p(±l/ (2A C ) ,to) = and p(0,to) = 
8A 2 , one can define the spatial width of the hump 
part in rogue wave as 1/A C . At the fixed time t = 

to, we have integral f^ 2 ^^ p(x,to)dx = AA C and 

/-i /(2Ac) P( x ^o)dx + fy(2A c ) P( x ^o)dx = -4A C . These 
results demonstrates clearly that for the attractive in- 
teratomic interaction, the generation of rogue wave with 
stronger breather characteristic is mainly due to the ac- 
cumulation of energy and atoms toward to its central 
part. The time- independent area relation shown in Fig. 
3(a), i.e., Si + 62 = S3, means that the loss of atoms 
in background completely transfer to the hump part of 
rogue wave. The forthcoming fundamental problem is 
that how rogue wave gather atoms and energy toward to 
its central part from the background. To this purpose, 
we investigate the exact number of atomic exchange be- 
tween rogue wave and background which has the form 



\ipR (x,t) - IpR (±00, t) I dx 



4ttA c 



y/l + AA A C T 2 ' 
(5) 

From the above expression, we can see that is time- 
aperiodic which is different from periodic exchange of 
atoms between the bright soliton and the cw background 
in Eq. ((2]). As shown in Fig. 3(b), the atoms in back- 
ground is gathered to the central part when t < to, which 
leads to the generation of a hump with two fillisters on 
the background along the space direction. The maxi- 
mal peak of the hump and the deepest fillisters occur at 
t = to- However, the atoms in the hump start to spread 
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FIG. 4. (Color online) (a) Evolution of the exact rogue wave 
solution Eq. (3]), where A c = 1, k c = 0, A = 0.01 and t = 10. 

(b) Evolution of the numerical solution of rogue wave with the 
initial condition Eq. ©, where A s = 1.92, A = 0.01, to = 10. 

(c) and (d) , The comparison of intensity profile between near- 
ideal rogue wave (in (c) with solid blue line) or sub-rogue wave 
pair (in (d) with solid blue line) and ideal rogue wave (circles). 



into the fillisters when t > to. Therefore, the hump grad- 
ually decay which verifies that the rogue wave is only one 
oscillation in time and displays a unstable dynamical be- 
havior. 

Considering the dynamics of rogue wave in the back- 
ground, on the one hand, a necessary condition for real- 
izing rogue wave in experiment is that the scale of rogue 
wave must be very small compared with the length of 
the background of BECs. In the real experiments [8[, the 
length of the background of BECs is at least 370/im. At 
the same time, in Fig. 3(a), the actual width of rogue 
wave is about 10a^ = 14/xra <C 370/im (a unity of coordi- 
nate corresponds to a± = ^h/muj± = 1.4/xm). Thus the 
rogue wave is observable experimentally. On the other 
hand, the decay rate of atoms in rogue wave can be con- 
trolled effectively by modulating the trapping frequency 
of external potential shown in Fig. 3(b). The decay time 
t d is about 8.05ms for A = 0.01 (uj± = 2tt x 700Hz 
and ujo — 2in x 7Hz originate from the experiment |9[), 
while t d ~ 2.30ms for A = 0.2 (u± = 2tt x 700Hz and 
ujo = Aiir x 70Hz), which demonstrate that a trap with 
small trapping frequency is conducive to the observation 
of matter rogue wave in BECs experiments. 

A question arises about the possibility of the creation 
of such a matter rogue wave experimentally. Generally, 
the excitation of rogue wave can be recovered by means 
of the numerical simulation with some particular initial 
conditions [29|. From the experimental points of view, 
optimal initial conditions are not only conducive to the 
experimental preparation, but it is also useful for under- 
standing the necessary conditions and physical mecha- 
nism of the generation of rogue wave. In what follows, 
we will look for the optimal initial conditions which can 
excite the resemble physical behavior of rogue wave. 



By comparing the Fig. 1(c) and (f) with Fig. 1(g), 
Eq. (|2]) with infinity oscillation period is a very good 
approximation of rogue wave solution Eq. (J3j), this pa- 
rameter regime yields characteristic rogue wave features 
in the spatio-temporal envelope even though the ideal 
rogue wave solution exists only asymptotically in the lim- 
itation of A c = 2A S and k c = k s . Based on modulation 
instability mechanism in ultracold atom system [30[, we 
consider the case of A 2 C > A 2 /A with k s = k c — cor- 
responding to the modulation instability process of cw 
background. In this case, we can take T « —to at t = 
for a very small quantity A, and the suitable values of 
A s and A c can ensure 6 « toA s M R /2 to be so large that 
k = e~ e is a small quantity. By linearizing the initial 
value with the small quantity n in Eq. ([2|) we get 

^(a?,0) = (a + 7 ^cosa )e^ (:c ' 0) , (6) 

where a = (2A\ - A 2 S - iA a M R )/(2A c ),<y = 
A s M R (M R -iA s ) /(2A 2 C ) and a = M R e~ xt °x, <p c « 
—A 2 to and the modulation frequency of initial condition 
Eq. (|6j) is Q = M R e~ xt °. The solution of initial value 
problem associated with Eq. (pQ) with initial condition 
Eq. J6l) can be described well by the exact solution Eq. 
([2]) [27]. For the creation of rogue wave, we require the 
value of A s /2 to approach A c associated with a very small 
modulation frequency in Eq. ©. The numerical results 
are shown in Fig. 4(b), which demonstrates that a small 
periodic perturbation with a very small modulation fre- 
quency can induce a near-ideal rogue wave localization, 
whose profile is basically consistent with the ideal the- 
oretical limit solution Eq. (|3]) as shown in Fig. 4(a) 
and (c). However, the obvious difference is that owing 
to the actions of modulation instability and instability of 
rogue wave, the initial near-ideal rogue wave can break up 
into two lower amplitude but equally strongly localized 
sub-rogue wave, and each sub-rogue wave itself exhibits 
ideal rogue wave characteristics as shown in Fig. 4(d), 
which agrees well with the optical experimental conclu- 
sions [2l|, [22[. As a result, a small initial periodic pertur- 
bation with a small modulation frequency can induce the 
generation of a near-ideal rogue wave by the modulation 
instability mechanism in BECs. 

Inspired by the experiments 0-[9(| , we can design the 
following experimental steps to observe the interesting 
rogue wave phenomenon in BECs: (i) Creating a conden- 
sates of 7 Li with total number of particles N ^ x 10 3 and 
continuous wave phase distribution by using quantum 
phase imprinting technique and amplitude engineering; 
(ii) Loading the condensates into a slightly expulsive har- 
monic potential with the parameters u± = 2tt x 700Hz 
and uo = 2in x 7Hz, and ramping up the scattering 
length in the form of a s (t) — — 0.9aBe At ; (hi) Optimal 
initial state Eq. (J6j) can be produced by imprinting a 
periodic perturbation laser with very small modulating 
frequency onto condensates. The main effect of this ex- 
pulsive term is that the center of the condensates accel- 
erates along the longitudinal direction. In addition, a 
crucial question is that since we require the scattering 
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length to change over time in above experimental proto- 
col, we must ensure the validity of quasi-one-dimensional 
regime and avoid the collapse of condensates with at- 
tractive interaction. In other words, we must ensure 
that the energy of two body interactions is much less 
than the kinetic energy in the transverse direction, i.e., 
e 2 ~ N\a s \/ao <C 1. With initial scattering length 
a s {t = 0) = -0.9a#, we can obtain e 2 « 0.03 <C 1. Af- 
ter about 50 dimensionless units of time, the scattering 
length becomes |a s (t)| = 1.4a# corresponds to e 2 ~ 0.057 
<C 1. So the validity of quasi-one-dimensional system can 
be maintained well. Finally, we emphasize that the inter- 
esting phenomenon of rogue wave can be observed within 
current experimental capability. 



IV. CONCLUSIONS 

In conclusion, we have investigated the formative 
mechanism of matter rogue wave in BECs with time- 
dependent interaction in an expulsive parabolic potential, 



analytically and numerically. The generation of rogue 
wave with stronger breathing characteristic is mainly due 
to the accumulation of energy and atoms toward to its 
central part. Rogue wave can not keeps dynamic sta- 
bility because of the aperiodic exchange of energy and 
atoms with the background; The decay rate of the num- 
ber of atoms in rogue wave can be controlled effectively 
by modulating the trapping frequency of the external po- 
tential. Our numerical results show that a small periodic 
perturbation with a smaller modulating frequency can in- 
duce the generation of the near-ideal rogue wave. Finally, 
we emphasize that the interesting phenomenon of rogue 
wave in BECs can be observed experimentally. Our work 
may facilitate the deeper studies of hydrodynamic rogue 
wave. 
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